Thermal spin-transfer in Fe-MgO-Fe tunnel junctions 
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We compute thermal spin transfer torques (TST) in Fe-MgO-Fe tunnel junctions using a first 
principles wave function-matching method. At room temperature, the TST in a junction with 3 
MgO monolayers amounts to 10 -7 J/m 2 / K, which is estimated to cause magnetization reversal for 
temperature differences over the barrier of the order of 10 K. The large TST can be explained by 
multiple scattering between interface states through ultrathin barriers. The angular dependence of 
the TST can be very skewed, possibly leading to thermally induced high-frequency generation. 

PACS numbers: 72.25.Ba, 85.75.-d, 72.10.Bg 



Spin-dependent thermoelectric effects in metallic mag- 
netic systems have been known for quite some time 
[l| but recently experience renewed interest. Spin 
caloritronic phenomena Q include the spin Seebeck ef- 
fect which should be distinguished from the spin- 
dependent Seebeck effect in nanostructures [4j|. Large 
spin-related Peltier cooling effects have been measured 
in magnetic NiCu nanopillars Hatami et al !6] pre- 
dicted that a temperature gradient induces a spin trans- 
fer torque that can excite a magnetization. Experimental 
evidence for the thermal spin-transfer torque has been 
obtained for Co-Cu-Co nanowires [7j. Slonczewski re- 
cently argued that thermal torques can be generated effi- 
ciently in spin valves with polarizing magnetic insulators 

1. 

Magnetic tunnel junctions (MTJs) of transition met- 
als with MgO barriers 0, [To| have great potential for ap- 
plications in magnetic random access memory (MRAM) 



Left lead Scattering region Right lead 




FIG. 1: Schematic Fe-MgO(3ML)-Fe(001) MTJs. We con- 
sider both a temperature difference AT and voltage difference 
Vb between the ferromagnetic reservoirs. The magnetization 
Mi of the left lead is fixed along the z-axis, while the mag- 
netization M2 of the right lead is rotated by an angle in 
the xz plane relative to Mi. The small dark gray (red and 
dark blue) squares in the scattering region represent O and 
Mg atoms, respectively, while the light gray (light blue) ones 
denote oxygen vacancies. 
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An im- 



elements and high-frequency generators 
portant goal of academic and corporate research remains 
the reduction of the critical currents necessary to induce 
magnetization precession and reversal 



dependent Seebeck effects in MTJs have very recently 
been computed [ItJ and measured [l8[ . A spin accumula- 
tions has been injected thermally into silicon by permal- 
loy contacts through MgO tunnel junctions [19j. 

Here we predict very large thermal spin transfer 
torques in MTJs with thin MgO barriers, which might 
open new possibilities to design memory elements and 
high-frequency generators driven by heat currents only. 
We have been motivated by the strong energy depen- 
dence of electron transmission through MTJs with thin 
barriers due to the existence of interface resonant states 
[20| , which should cause large thermoelectric effects. Fo- 
cussing on epitaxial Fe-MgO-Fe MTJs under a temper- 
ature bias, we demonstrate the effectiveness of thermal 
spin transfer torques by ah initio calculations based on 
the Landauer-Biittiker transport formalism. Consider an 



MTJ as sketched in Fig. [I] with a voltage and temper- 
ature bias over the two leads, which are in local ther- 
mal equilibrium with Fermi-Dirac distribution functions 
IS [16|. Spin- f L/R ( e ) = [ e ^-^/R)/k B T L/R + l]- 1 and local chemical 
potentials Hl/r an d temperatures T L j R . The generalized 
Landauer-Biittiker formalism [2l[ is very suitable to han- 
dle transport through layered magnetic structures. The 
spin current from the n-th layer to the n + 1-th layer can 
be written as 



22] 



Jn+i,n = j de [t£ +ljn (e)/ L (e) + t* n+1 (e)f R (e)] . 

(1) 

Here the energy-dependent spin transmission coeffi- 
cient matrix from the left (right) direction is defined 



L/R 
L n+l,n 



E kll (* 



L/R, 



All \«k„ 'tofci,»(k||)|<*(e)> 



with spin current operator 7U7 n +i,n(k||) = 
_ Re {^#nL,n+ii/ (k|| ) j , H n L,n+iL' (k|| ) de- 

notes the Hamiltonian matrix in spin space [22|, where 
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L = (Z,ra) are the azimuthal and magnetic quantum 
numbers, a is the vector of Pauli spin matrices and ky 
is integrated over the two-dimensional Brillouin zone of 
transverse modes. 

When the applied voltage vanishes and A^AT < ej? 
we may expand Eq. $Q in AT = Tr — Tl- 

J n+i,n = ^{ J ^e/(eF,To)[t^ +1?n (e) + t£ n+1 (e)] 

+ Wo I de{€ ~ e 4 [t ' +1 ' n(e) " t -'-+ l(e)]} ' 

(2) 

where T = (T L + T fl )/2 . The first term in J n+ i, n 
= Jn+i n + Jn+i,n * s the equilibrium spin current that 
communicates the exchange coupling through the barrier, 
while Jn+i,n i s the thermal spin current. The torque 
acting on the n-th layer is the difference between the 
incoming and outgoing spin currents T^ T = J^-i ~ 
Jn+i n- The total thermal spin transfer (TST) torque is 
then 



Ta7 



AT 

Wo 



d 

de(e - e F )rv(e) — f(e), (3) 



where r v (e) = (e/167r) E^Lw^n-ifc) ~ t n-i,n( e ) - 
^n+i n( e ) + n+i( e )] * s the electrical torkance, and the 
sum from N to oo runs from the first layer at the in- 
terface until deep into the bulk of the "free" magnetic 
lead. The thermal torkance tt = Tat/ AT depends on 
the energy dependent transmission and To. Only when 
Ty{e) varies slowly around the Fermi level in the ther- 
mal window fc^AT, the Sommerfeld expansion may be 
employed and r T — >■ — (e/2) (7r 2 /c^T /6) dry (e)\ 6F . For 
comparison, the linear-response volt age- driven torkance 
reads [22[ T av (ep) /Vb — >■ Tv(ei?), where is the ap- 
plied bias. We found that Eq. (j3j) accurately reproduces 
non-linear calculations based on Eq. (pQ) in the parameter 
regime considered here. 

Here we focus mainly on Fe-MgO-Fe(lOO) with 3 lay- 
ers (3T) MgO (<~ 6 A), corresponding to the thinnest 
barrier that can be reliably grown gnorins; the 

small lattice mismatch between the lead and barriers, we 
assume an isomorphous structure with unrelaxed interfa- 
cial atoms at bcc positions. We use a 1200 x 1200 fc-point 
mesh in the full two-dimensional Brillouin zone (BZ) of 
transport channels to ensure numerical convergence. De- 
tails of the electronic structure and transport calculations 
can be found elsewhere [23| . 

The energy-dependent transmission through our MTJ 
is plotted in Fig. [2] for the parallel (P) and antiparal- 
lel (AP) co nfig urations. Vacancies break the crystalline 
symmetry [24[, broaden the resonant peaks of the mi- 
nority spin channel, and enhance the coupling between 
majority Ai state of one lead and the surface states of 
the minority spin on the other side. The newly opened 
channels increase the AP conductance. 
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FIG. 2: Energy-dependent transmission of Fe-MgO(3L)- 
Fe(001) MTJs with (a) perfect interfaces and (b) 10% OVs 
at both interfaces. The shaded area indicates the thermal 
energy window ksTo at room temperature. 



For specular Fe-MgO interfaces, we find a zero- 
bias "optimistic" tunneling magnetoresistance (TMR) of 
1300% and a resistance area 0.063 Q jim 2 . Our TMR ra- 
tio is consistent with a "pessimistic" TMR calculated to 
be around 0.93 for the same barrier thickness [25| with 
majority-spin transmission around 0Ae 2 /h jifj. When 
10% oxygen vacancies (OVs) (the energetically most 
favorable defects) are introduced at the interface, the 
TMR decreases to 96% and RA = 0.036 Qfim 2 . When 
10% OVs are introduced in the middle of barrier, the 
TMR drops to 70%, and the RA slightly increases to 
0.076 Q /im 2 . These results are comparable with the mea- 
sured 0.19 Q /im 2 . and TMR = 15% for a similar barrier 
thickness at room temperature [l3[ . 

From the energy-dependent transmissions, we can also 
compute the magneto Seebeck coefficients [17] and elec- 
tronic heat conductances. At T = 10 K, we find = 
2.1 x 10 6 WK _1 m- 2 and = 0.14 x 10 6 WK _1 m- 2 , 
which, including estimated phonon contributions to the 
heat conductance, leads to the thermoelectric figure of 



merit (ZT) 10K ~ lO" 3 . 

In Fig. [3] (a) and (b), we present the in- and 
out-of-plane angular-resolved torkances of specular Fe- 
MgO(3L)-Fe MTJs. The in-plane torkance is smooth 
in most energy regions, indicating good numerical con- 
vergence. We observe two resonances: a small one at 
E F -0.02 eV that contributes to the TST for T > 100 K 
and a sharp and larger peak at Ep — 0.0725 eV that 
contributes to the TST for T > 300 K. The out-of- 
plane torkance is much more sensitive to numerical de- 
tails. The noise in Fig. [3] does not affect the integrated 
TST, however. At energies far away from the Fermi level 
(E > E F + 0.03 eV and E < E F - 0.09 eV), the in-plane 
torkance is small and proportional to sin# as predicted 
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FIG. 3: Energy and angle-dependent (a) in-plane and (b) 
out-of-plane torkance, angular dependent (c) in-plane and (d) 
out-of-planeTST, and ratio of (e) in-plane and (f) out-of-plane 
TST to thermocurrent (Tat /I at) of epitaxial Fe-MgO(3L)- 
Fe(OOl) MTJs. 

by model studies [27]. However, this region contributes 
only weakly to the TST. The two sharp peaks near the 
Fermi level show an angular dependence that deviates 
strongly from a sine function. The asymmetry of the 
angular dependence of the in-plane TST reflects multi- 
ple scattering in the barrier and is therefore exponentially 
suppressed for thick layers. The in-plane TST of 7L MTJ 
(not shown) already agrees well with a sine function. 

The angular dependence of the observable TST, i.e. 
the energy integral in Eq. (j3j), is plotted in Figs. [3^c) and 
(d) . We observe strong deviations from a sine function at 
all temperatures considered. The skewness can be traced 
to multiple-reflection hot-spots caused by the interfacial 
resonances mentioned above. At room temperature the 
in-plane torkance peaks around 165° and the functional 
form can be fitted to an asymmetry parameter [28| of 
A = 3.5. This value is much larger than observed for 
the voltage-induced torque in metallic spin valves [29| . 
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FIG. 4: k|| -resolved electron transmission probability of (a) 
majority- and (b) minority-spin, and in-plane torkance from 
(c) majority- and (d) minority-spin channels in 3L MTJs with 
magnetization angle of 177.5° at energy Ef — 0.0725 eV. Ma- 
jority (minority) are defined for the left lead. The integrated 
transmission probabilities are 4.6(6.3) x 10 -3 e 2 /h and inte- 
grated in-plane torkances are 5.5 x 10 ~ 5 (—4.2 x 10 _3 )eV for 
majority (minority) spins, respectively. 

which should be beneficial for high-frequency generation 
. We therefore suggest the possibility of efficient spin 
oscillators driven by heat flows through MTJs. The out- 
of-plane term is an effective field that dominates the in- 
plane term for angles > 165°. 

Fig. G3e) and (f) displays the angular dependence of 
the spin transfer efficiency monitored by the ratio of the 
torque to charge current density for a given temperature 
bias AT = 1 K. We find that the ratio (for both in-plane 
and out-of-plane terms) increases strongly close to the 
APC, for which the charge current is suppressed by the 
high spin polarization of the Fe-MgO interface (30| . 

The high spin transfer efficiency near the APC can be 
explained by multiple reflection due to resonant tunnel- 
ing. In Fig. [4] resonant tunneling is observed in the APC 
at a chosen energy in both spin channels with a con- 
ductance polarization of 16%. Their contributions to the 
torkance are much larger, since the minority spin channel 
transfers 99% of the torkance due to its high interfacial 
electronic density of states. Here majority and minority 
spins are defined for the left lead. The resonance persists 
in the exact APC (Fig. Efa)) but spin transfer vanishes 
for collinear magnetizations. 

In Table U we compare TSTs equivalent to AT = 1 K 
at To = 300 K with electric STs for MTJs with 3, 5, and 
7 MgO layers. The equivalent bias and current density 
of the thinnest barrier sample is much larger than that of 
the thicker one, which reflects the exponential decay of 
the conductance as a function of barrier thickness. AV eq 



TABLE I: Thermal torque Tik per unit cell in nL MgO 
MTJs at T = 300 K and AT = 1 K under closed and 
open (in brackets) circuit conditions for = 90°. A14 q = 
Tik/t v is the equivalent bias. AVik/Ak is the thermovolt- 
age/thermocurrent, TV and iy are electrically induced torque 
and current, respectively. 

AV eq AV 1K T 1k /Iik Tv/Iv 



(mJ/V/m 2 ) 
3 0.72 
5 0.082 
7 0.011 



Tik* 
(nJ 

-195(-232) -0.27 
-3.32(-5.33) -0.040 
-0.24(-0.062) -0.021 



IK AVeq AVlK ± IK / llK '±V / lv 

im" 2 ) (mV) (mV) (h/2e) (h/2e) 
\(-9.Z9.\ -0 97 059 -0 Q4 91 



0.052 -0.94 
0.025 -0.84 
-0.0154 3.58 



0.21 
0.46 
0.98 



n jy-'m-' 2 = 3 x 10^/2) kn _1 m 



is the ratio of thermal to electric torkance, which is larger 
for 3-MgO (IK - 0.27 meV). AV 1K demonstrates that 
TSTs decrease faster than the electric STs when the bar- 
rier gets thicker. The sign change in AVik as a function 
of barrier thickness is attributed to that of the Seebeck 
coefficient. Moreover, the torque to current density ratio 
T/I is larger for the thermal than the electric case, in- 
dicating the superior efficiency of spin angular moment 
transfer by temperature differences. 

The TST is potentially useful for manipulating the 
magnetic configurations in MTJs with thin barriers. We 
estimate the critical temperature bias AT C by comparing 
the TST with the measured torques at the critical volt- 
age biases in CoFeB MTJ's at room temperature (3l| . 
For 3ML MgO the switch from APC to PC should oc- 
cur close to AT C AP ^ P = 6.5 K since then |T AT /sin<9| = 
20 x 10 -6 Jm~ 2 equals the critical torque for electric 
switching. At AT C P ^ AP = 56.5 K, |T AT /sin(9| = 8.2 x 
10 -6 Jm~ 2 equals the critical torque for electric PC to 
APC switching (3l|. We note that tt is function of the 
global temperature that saturates around 275 K. Room 
temperature conditions are therefore favorable for ther- 
mal magnetization switching. 

In an open circuit, the thermoelectric current vanishes, 
but not the thermospin current, thereby allowing transfer 
of angular momentum without transfer of charge. The 
thermal torque is even found to be larger in the closed 
compared to the open circuit, since the equivalent bias 
AV e(l and the thermovoltage AVik have opposite signs. 

The spectral features due to resonances are sensitive 
to disorder. In Fig. [5] we show the angular-dependent 
torkance in 3-MgO with 10% OVs at both interfaces. We 
make comparison for two situations: one is at Ef, and 
another is at resonant peaks near to Ef- The resonant 
peaks in the clean samples at Ep — 0.0725 eV shift to 
lower energy (around Ep — 0.055 eV) in the presence of 
10% OV as shown in Fig. [2j so different energies are cho- 
sen to compare clean and dirty situations. We observe 
that the disorder to a large extent restores the sin 6 angu- 
lar dependence. The order of magnitude of the in-plane 
torkance of the ideal junctions at Ep is unmodified. The 
situation at the resonant peak is more complicated by a 
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FIG. 5: Angular dependent (a) in-plane and (b) out-of-plane 
torkance in Fe-MgO(3L)-Fe(001) MTJs with specular and dis- 
ordered interfaces at two selected energies, i.e., at the Fermi 
energy and at the resonance. Squares and circles are results 
for specular interfaces at Ef and Ef — 0.0725 eV, respec- 
tively; the up-pointing triangles and down-pointing triangles 
are disordered samples at Ef and i^F — 0.055 eV, respectively. 



shift from E F - 0.0725 eV (specular) to E F - 0.055 eV 
(disordered) with decreased amplitude. A full calculation 
of the TST in the presence of OV disorder at room tem- 
perature is beyond our present computational capacity, 
but the two noted changes of the resonance will at least 
partly cancel each other. 

In summary, we calculate TSTs of the order of 
10 -6 Jm" 2 in ultrathin Fe-MgO-Fe tunnel junctions at 
room temperature for a temperature bias of 10 K. A 
strong asymmetric angular dependence of TSTs is pre- 
dicted for ballistic junctions. Based on these results 
we predict heat-flow-induced magnetization reversal and 
high frequency generation in magnetic tunnel junctions. 
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